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Abstract 

We study the optimisation of exact renormalisation group (ERG) flows. We explain 
why the convergence of approximate solutions towards the physical theory is opti- 
mised by appropriate choices of the regularisation. We consider specific optimised 
regulators for bosonic and fermionic fields and compare the optimised ERG flows 
with generic ones. This is done up to second order in the derivative expansion at 
both vanishing and non-vanishing temperature. We find that optimised flows at 
finite temperature factorise. This corresponds to the disentangling of thermal and 
quantum fluctuations. A similar factorisation is found at second order in the deriva- 
tive expansion. The corresponding optimised flow for a "proper-time renormalisation 
group" is also provided to leading order in the derivative expansion. 
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I. INTRODUCTION 



Wilsonian renormalisation group techniques |l],|2| like the exact renormalisation group 
(ERG) are important tools for addressing non-perturbative problems within quantum 
field theory (for recent reviews, see Refs. [|7|f|). They are similar in spirit to the block-spin 
action invented in condensed matter physics, and their particular strength is their flexibility, 
allowing for systematic approximations without being tied to the small coupling region. The 
ERG is based on an infrared (IR) regularisation with momentum scale parameter k of the 
full propagator, which turns the corresponding effective action into a scale dependent func- 
tional Tfc. The ERG flow describes the change of the effective action under an infinitesimal 
variation of the IR scale k. It thereby interpolates between the initial UV action r^ = A and 
the full quantum effective action T = T^q. Although the flow depends explicitly on the 
specific infrared regulator chosen, the endpoint of the integrated full flow does not. 

An explicit computation of the IR effective theory based on the ERG flow requires the 
specification of the field content, the initial condition I\ and the choice of a particular IR 
regulator. The UV initial condition is typically given by the classical action. Hence, the 
main physical information is contained in the ERG flow itself. Most problems of physical 
interest are too complex to be solved exactly and an application of this formalism - as of 
any other method - is bound to certain approximations. Furthermore, the flow equation is 
equivalent to infinitely many coupled partial differential equations, which would seem very 
difficult to solve exactly. Therefore, one has to resort to some approximations or truncations 
which allow, at least in principle, for a systematic computation of the full quantum effective 
action. In order to provide reliable physical predictions, like a high precision computation 
of universal critical exponents, it is mandatory to provide a good control for approximated 
ERG flows. 

A number of systematic expansion schemes for flow equations are known, including stan- 
dard perturbation theory. Non-perturbative expansions of the effective action, not bound to 
the weakly coupled regime, are the derivative expansion, expansions in powers of the fields, 
or combinations thereof. For example, the leading order of the derivative expansion retains 
only an effective potential and a standard kinetic term, and contains non-perturbative infor- 
mation as it corresponds to the resummation of infinitely many perturbative loop diagrams. 
The study of approximate quantum effective actions along these lines is a sensible procedure 
since the underlying expansions admit a systematic improvement to higher order. 

Solutions to truncated flow equations display a spurious dependence on the IR regulator 
P^n|. This is similar to the scheme dependence of physical observables observed within 
perturbative QCD [ljj, or the truncation dependence of solutions to Schwinger- Dyson equa- 
tions. Its origin is the following. The IR regulator couples, through the flow equation, to 
all vertex functions of the theory. The flow trajectory of the functional in the space 
of all effective action functionals depends on the regulator. Hence, the regulator - while 
regulating the flow - also modifies the effective interactions at intermediate scales k ^ 0. 
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In other words, the effective action at intermediate scales still has some memory of the 
details of how the integrating-out of degrees of freedom has been performed. This regulator 
dependence is of no relevance for the full flow. Eventually, the convergence towards the full 
quantum effective action for any regulator ensures that all regulator-induced interactions 
cancel out in the physical limit. Approximations imply that certain vertex functions are 
neglected. Then, not all regulator-induced interactions cancel out for k — > 0: the missing 
back coupling of the neglected vertex functions is responsible for regulator-dependent terms 
in the physical limit. In consequence, approximations to the full quantum effective action 
depend spuriously on the scheme. 

Recently, a new line of reasoning has been put forward which essentially turns this obser- 
vation around [IS]: given that the solution of a truncated flow depends on the regulator, it 



should be possible to identify specific ones which "optimise" the physical content of a given 
approximation. Optimised regulators stabilise the flow and lead to a faster convergence of 
expansions, such that the main physical information is almost exclusively contained within 
a few leading terms, and higher order contributions remain small [18] . 



In Ref. [|L5 |, we have derived a simple and generic optimisation criterion for ERG flows, 



based only on the full inverse propagator at vanishing field. Given the set of possible IR 
regulators, the criterion allows to distinguish the quality of regulators in the sense outlined 
above. In the present paper we study a specific "optimised" regulator for both bosonic 
and fermionic degrees of freedom. To be more explicit, we introduce the ERG flow for the 
effective action [§-0. For bosonic fields 0, it is given by 

9 ' r ^=^( wShj +Rt ) (L1) 

Here, the trace denotes a sum over all loop momenta and indices, and t = In A; is the 
logarithmic scale parameter. The flow has a simple one-loop structure. The Wilsonian 
"integrating-out" is achieved by the infrared regulator Rk- It regulates the propagator for 
small momenta, while the insertion d t Rk cuts off the large-momentum contributions. In 
total, only a small momentum window about q 2 m k 2 contributes to the flow. Apart from a 
few constraints displayed later, the function can be chosen at will. A "good" choice for 
the regulator function is at the root for reliable physical predictions, and we consider, for 
the bosonic fields, the optimised regulator 

RT(q 2 ) = Z k -(k 2 -q 2 )e(k 2 -q 2 ), (1.2) 

where is an appropriately defined wave function renormalisation. This regulator is par- 
ticularly simple: for loop momenta q 2 > k 2 it vanishes identically and the effective prop- 
agator appearing in the flow equation is not modified; for loop-momenta q 2 < k 2 it acts 
like a momentum-dependent mass term in such a way that the inverse effective propagator 
~ q 2 + Rk(q 2 ) becomes a momentum independent constant. In consequence, the effective 
infra-red propagator does no longer distinguish between the different modes with q 2 < k 2 . 
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Optimised flows based on Eq. (|1.2| ) derive from a generic optimisation criterion | 15| , and 



have a number of remarkable properties. The optimised flow leads to the fastest decou- 



pling of heavy modes, in accordance with the decoupling theorem [19|. In the limit k — > 



optimised flows smoothly approach a convex effective action, owing to a simple analytic 
pole of the flow [|15|,|]1|. At non- vanishing temperature, the optimised flow factorises: the 
contributions from thermal and quantum fluctuations are disentangled, unlike for generic 
flows. A similar factorisation of the flow holds to second order in the derivative expansion 
for field-independent wave function renormalisations, and a partial factorisation is found 
for the general case. Finally, the optimised flow has a very simple analytic structure. This 
facilitates their study and is helpful for both analytical or numerical considerations. All 
these properties lead to a stabilisation of the flow and an improved convergence towards the 
physical theory. Analogous results for fermionic flows are discussed as well. 

We also study this question within an RG formalism based on a proper-time regularisa- 
tion of the operator trace for the one-loop effective action [2Q|, which we call "proper-time 
renormalisation group" (PTRG) for short. In contrast to the ERG, the PTRG has no path 



integral derivation, which makes the conceptual reasoning more difficult [21]. Still, owing 



to the close similarity to the ERG at leading order in the derivative expansion, it is possible 
to identify the analogue of Eq. (|1.2|) for the PTRG. 

The format of the paper is as follows. We introduce the physical ideas behind the generic 
optimisation condition. Explicit realisations for bosonic and fermionic degrees of freedom are 
introduced as well (Section [TT|). The main characteristics of optimised flows are discussed 
to leading order in the derivative expansion, and contrasted with those of generic flows 
(Section PTI]). We then turn to the discussion of quantum field theories at finite temperature. 
We show that optimised thermal flows factorise on the level of the flow equation, unlike 
generic flows. A simple physical explanation for the factorisation is provided (Section [TV| ). 
Next, we consider the extension to higher orders in the derivative expansion. The cases 
of field dependent or independent wave-function renormalisations are both discussed, and 
a similar factorisation for optimised flows is established (Section |V|). Finally, we provide 
the corresponding optimised proper-time cut-off for the PTRG (Section [Vf). Due to the 
qualitative difference of the topics studied, we discuss our findings separately at the end of 
the corresponding sections. We close with a summary and an outlook (Section |VII|) . Three 
Appendices contain technical details and explicit expressions for optimised flows. 



II. OPTIMISATION 

In this section, we discuss a generic optimisation criterion for ERG flows for Euclidean 
quantum field theories. In particular, we provide a simple and explicit optimised regulator 
for both bosonic and fermionic flows. Prior to this, we have to review a few basic properties 
of IR regulator functions, which are at the root of the subsequent considerations. 
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A. Regulators 



The flow equation ( |1.1|) is denned through the infrared regulator functions Rk(q 2 ) and 
RF,k(q 2 ), respectively These operators depend on an infrared scale k, which induces 

a scale dependence. When written in terms of the scale-dependent effective action Tk, 



the scale dependence is given precisely by the flow equation (IT). The right-hand side of 
Eq. (|1T|) contains the full inverse propagators and the trace denotes a sum over all indices 
and integration over all momenta. 

The regulator scheme (RS) functions can be chosen at will, however, within some basic 
restrictions. These restrictions ensure that the flow equation is well-defined, thereby inter- 
polating between an initial action in the UV and the full quantum effective action in the 
IR. More specifically, it is required that 

lim R k (q 2 )>0. (2.1) 

This ensures that the effective propagator at vanishing field remains finite in the infrared 
limit q 2 — > 0, and no infrared divergences are encountered in the presence of massless 
modes. This property makes Rk an infrared regulator. If the limit ( |2.1|) is finite, we call 
the corresponding regulator masslike. The second requirement is the vanishing of Rk in the 
infrared, 

lim R k (q 2 )^0. (2.2) 

This guarantees that the regulator function is removed in the physical limit, where the 
scale-dependent effective action T k reduces to the quantum effective action T = lim fc ^ r fc . 
The third condition to be met is 

lim^g 2 ) -> oo . (2.3) 

This way it is ensured that Tk approaches the microscopic action S = lim.k^,A^k in the UV 
limit k — * A. In the rest of the paper, we put A = oo for the UV scale, although our main 
line of reasoning can be applied for finite A as well. With this choice, the regulator function 
depends only on q 2 and k 2 , and it is convenient to introduce a dimensionless function 

r(q 2 /k 2 ) as 

R k {q 2 ) = Z k q 2 r{q 2 /k 2 ) . (2.4) 

with Zk an appropriate wave function renormalisation (cf. Section Z k = 1 to leading 
order in the derivative expansion. Owing to the general conditions imposed on the regulator, 
the function r(y) ranges between < r(y) < oo. 

Another condition concerns the proper normalisation of the regulator. The normalisation 
fixes the scale at which the IR regulator becomes effective. Let us employ the condition 
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Rk{q 2 = c B k 2 ) = Z k c B k 2 



(2.5) 



for bosons (a similar condition holds for fermions, see below) and cb > 0.[] The normali- 
sation translates into the condition r(ce) = 1. Two different choices for cb can always be 
mapped onto each other through a rescaling of the IR scale k. Hence, a proper normal- 
isation is only of relevance for a comparison of different regulators (as done in Ref. ||15||), 
or for theories containing different bosonic and/or fermionic degrees of freedom, where the 
relative normalisation of the regulators can become important. 



B. Optimisation criterion 

Here, we discuss an optimisation criterion for ERG flows, which ensures that flows like 
Eq. ( |1 . 1|) and approximations to it have good convergence and stability properties. Follow- 
ing Ref. [jn| (see also Ref. [18|), we first provide the general criterion for optimised choices of 



RS functions. Then, more specifically, we apply this idea to bosonic and fermionic theories 
with standard kinetic terms. 

The physical information of the flow equation ( |1.1| ) is contained in the full effective 
inverse propagator, which is given by 

(2.6) 



5<j>(q)54>(-q) 



Notice that Eq. ( [2.6|) depends both on the fields and on the RS function. The ERG flow is 
well-defined as long as the full inverse propagator displays a gap, 



q2>6 y8(j){q)5(j)(-q) 



mm 




Ck 2 >0. (2.7) 



The functional derivative is evaluated at a properly chosen expansion point <j)Q. The exis- 
tence of the gap C > implies an IR regularisation. Furthermore, the gap is a prerequisite 



for the ERG formalism. Elsewise, Eq. (LI) becomes singular at points where the full inverse 
effective propagator develops zero modes|] The size of the gap C in Eq. ( |2.7| ) depends both 
on the RS function and on dimensionless parameters like (fl/k 2 or mass ratios, specific to 
the particular theory studied. 



tin Ref. [15 the convention cb = 1 has been used. 

*The case C = indicates that a saddle point expansion about 4>o is not applicable. Those points 
4>q in field space with C = correspond to an instability. The problem can be cured by chosing a 
more appropriate expansion point such that C > 0. For related literature, see Ref. |22|. 
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A natural optimisation criterion based on Eq. ( |2.7| ) consists of maximising the gap C 
over the space of all possible RS functions. Optimised RS functions are those for which the 
maximum of C is attained. The optimisation ensures that the momentum-dependent kernel 
of the ERG flow is the most regular. Therefore we expect that optimised flows are much 
more stable against approximations and show better convergence properties. 

The optimisation condition as formulated above is, essentially, only sensitive to the mo- 
mentum dependence of the full inverse propagator. Dropping momentum-independent terms 
on the left-hand side of Eq. (|2.7| ) changes the number C accordingly, but leaves the explicit 
dependence on Rk(q 2 ) unchanged. Therefore, the optimisation leads to the same set of opti- 
mised RS functions as long as the implicit dependence of rj^ [</>] = 5 2 Tk[(j))/5(f)(q)5(f)(— q) 
on the RS function remains negligible. For this reason, the optimisation condition of 



Refs. (T^JT^] is based only on the momentum- dependent terms of Eq. ( [2.6| ). 

From now on, we concentrate on a standard kinetic term. The effect of a field-dependent 
wave function renormalisation can be taken into account as well (see Section [V| below) . We 
expand the full inverse propagator as Zk{q 2 + Z^R^q 2 ) + . . .) about the regularised kinetic 



term. Finally, dropping the momentum-independent terms transforms Eq. ( |2.7|) into 

mm (q 2 + Z^ 1 R k (q 2 )) = C k 2 > . (2.8) 



A far reaching consequence of the infrared regulator in Eq. (|2.8|) is the presence of a gap for 
all k > 0, which follows trivially from Eq. ( j2.1|) . The decisive difference between Eq. ( |2.7| ) 
and Eq. ( |2.8|) is that the size of the gap C > in Eq. (|2.8| ) depends only on the particular 
choice for the RS, but not on the specific theory. Rewriting Eq. ( |2.8| ) in dimensions of k 
leads to 

P 2 {y) = q 2 /k 2 + R k {q 2 )/{Z k k 2 ) = y[l + r(y)} , (2.9) 
where y = q 2 /k 2 . Expressed in terms of Eq. (|2.9j), the size of the gap is given by 

C = mm P 2 (y) . (2.10) 

Any RS function is now characterised by the associated gap C. The size of the gap can 
be made arbitrarily small. Effectively, this corresponds to removing the IR regulator in the 
first place. However, for fixed normalisation Cb, it cannot be made arbitrarily large, C < oo. 
Hence, the natural optimisation condition, which is the requirement to maximise the gap, 
becomes 

^ Pt = max(minP%)). (2.11) 

A few comments are in order. The maximum in Eq. ( |2.11| ) is taken over the (infinite- 
dimensional) space of all possible RS functions. The number C opt is uniquely determined 
and reads C opt = 2c#, where c# is the normalisation of bosonic regulators. From now on, we 
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refer to Eq. ( [2. lip as an "optimisation condition" , and all RS functions for which C = C op t 



are called solutions to the optimisation condition. The space of solutions to the optimisation 
condition is infinite-dimensional. Notice also that the condition to minimise the gap is not 
an extremisation linked to the regulator, because it corresponds to removing the IR regu- 



larisation. In Ref . |T5| , a variety of different solutions have been found, and some examples 
are given in Fig. 1 below. 

In order to obtain Eq. ( |2.S| ), we have assumed a standard kinetic term for the fields. 
Therefore, the resulting optimisation condition Eq. ([2.1 1|) is independent on the specific 
theory. Once the momentum-dependent part of 1^ depends on the fields, the correspond- 
ing optimisation condition based on the momentum-dependent part of Eq. ( |2.7| ) is sensitive 
to the specific theory. Within a derivative expansion, this happens starting from the second 
order (cf. the discussion in Section |V|). 

The optimisation condition has a number of interpretations in more physical terms 



(cf. Refs. ||15| , |18ll ) . It has been shown that the radius of convergence for amplitude ex- 
pansions is given by C. Therefore the optimisation condition improves their convergence. 
Furthermore, it leads to a smooth approach towards a convex effective potential in the IR 
limit k —>■ 0. It has also been shown that it improves the convergence of the derivative 
expansion |18[ . Finally, it is worth emphasizing that the optimisation criterion is a rather 
mild condition: all regulator functions are described by at most countably infinitely many 
parameters, because Rk is at least square integrable. Of these, only one parameter is fixed 
by the optimisation criterion. 



We now turn to the discussion of fermionic degrees of freedom ijj and ip [f23 , 23 ] . The flow 
equation is given by 

WW, A = "Tr { M ^,; ( _ q) + 1 0*** ■ (2-12) 

As usual, the trace sums over all loop momenta and indices. The constraints on the function 
R F k are similar to those on ||24|| . Following Ref. |j24 |, we chose the regulator proportional 



to <f[ and introduce 

R F>k (q) = Z F)k i r F (q 2 /k 2 ) (2.13) 
We chose the normalisation as 

R 2 Fk (q 2 = c F k 2 ) = c F k 2 . (2.14) 

This translates into the condition r F (c F ) = 1. It has been shown that the fermionic analogue 
of the function Eq. ( |2.9| ) is given by |23| 



p 2 (y) = y[i+r F (y)] 2 . (2.15) 

Therefore, we can define the fermionic gap as 
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C F = min P 2 F (y), (2.16) 
y>0 

and the corresponding optimisation condition reads 

C F , opt = max ( min P F (y) ) . (2.17) 

(RS) \y>0 / 

The optimised fermionic gap is uniquely determined through the normalisation cp as Cp opt = 
4cp. Conceptually, the fermionic case is treated in the same way as the bosonic case. The 
sole difference stems from the fact that the bosonic kinetic term contains two derivatives, 
while the fermionic kinetic term contains only one. Therefore, the functions (|2.9|) and ( |2.15| ) 
entering the optimisation condition are different. 



C. Derivation of optimised bosonic and fermionic regulators 

A lot of effort has been spent in order to provide explicit regulators which lead to suf- 
ficiently simple and analytic ERG flows. For example, the sharp cut-off provides a simple 
explicit flow to leading order in the derivative expansion. For this reason, it is one of the 
most intensively studied flows in the field (cf. Refs. [p],p|,p5|,|6|, p6| , p7| ) . Other attempts have 



or 



been made based on power-like regulators Rk ~ q 2 (k 2 /q 2 ) b for b = 1 and b = 2 |28 
variants of a mass-term regulator Rk ~ k 2 Q(k 2 — q 2 ). These regulators are still sufficiently 
simple from an algebraic point of view, and lead to reasonably simple flows.0 However, in 
absence of an underlying "guiding principle" it was not obvious how to make progress given 
the plethora of possible regulators, and in particular, how to distinguish the "quality" of 
the corresponding flows. 

Here, in turn, we take full advantage of the existence of a guideline provided by the 
optimisation criterion. We propose a regulator which (i) solves the optimisation criterion, 
(ii) is based on an additional stability criterion for approximate flows, and (iii) leads to simple 
explicit expressions for the corresponding flows. The heuristic derivation runs as follows. 
The space of regulators which solve the optimisation criterion is still infinite dimensional. 
Let us seek for a "simple" solution to Eq. ( [2.11|) . The simplest one corresponds to an inverse 



propagator which is flat, i.e. momentum-independent, P 2 = C opt . Take C opt = 1. This 
immediately implies, using Eqs. (|2.4j ) and (2J3), that Rk{q 2 ) = k 2 — q 2 . Our naive Ansatz 



is consistent with Eqs. (|2.1|) and (|2.3|) , but not with the main requirement Eq. ( |2.2j ) for 
small k 2 < q 2 . In order to fulfil Eq. ( [2.2|) , the regulator has to be cut-off above some loop 
momenta. Therefore, a natural proposal for the bosonic case consists in taking 



RT(q 2 ) = (k 2 -q 2 )Q(k 2 -q 2 ) . (2.18) 



^Of these, only the power- like regulator with 6 = 2 solves the optimisation condition Eq. ( 2.11| ). 
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The ultraviolet modes q 2 > k 2 are not touched by this regulator because Eq. ( |2.18| ) van- 
ishes identically for q 2 > k 2 . In turn, for all modes with q 2 < k 2 the regulator acts as a 
momentum-dependent mass term ~ (k 2 — q 2 ) with the infrared limit ~ k 2 for vanishing mo- 
menta. It is a masslike regulator. By construction, the inverse propagator at vanishing field 
Eq. ( |2.6| ) becomes momentum independent for all q 2 < k 2 (see Fig. 1). It is this property 
which is responsible for the main characteristics of the regulator: all infrared momentum 
modes below the scale k are treated in the same way since the effective inverse propagator 
does no longer distinguish between them. 




0.5 



0.5 



1 

q 2 /k 2 



1.5 



Figure 1: Optimised inverse propagators P 2 pt for different regulators, normalised as 
r(^) = 1. The regulator Eq. ( 2.18| ) is given by the full line. The thin dashed line corre- 
sponds to r = 0. All other dashed lines, given for comparison, correspond to the different 
optimised regulators of Fig. 3 in Ref. [15|. 



The dimensionless regulator function r(y) is defined in Eq. (|2.4|) . With the choice made 
in Eq. (|2~T^ ) it follows that 

r op t(y)= ^-l)e(l-y). (2.19) 

The regulator function is normalised with cb = \- Such a normalisation can always be 
achieved. Other normalisations are considered at the end of Section |TTT| . 

In Fig. 1, we have displayed the effective inverse propagator at vanishing field for dif- 
ferent optimised regulators. The full line represents Eq. ( |2.18| ), and the thin dashed line 
corresponds to R = (no regulator). The set of dashed lines correspond to the optimised 
regulators discussed in Fig. 3 of Ref. [115] . Here, they have been given for comparison. Notice 
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that all curves cross in the normalisation point r(ce) = 1. All optimised propagators display 
the same gap C opt = 2cb, but differ essentially in the curvature around their minima. 

The fermionic analogue of Eq. ( |2.18| ) is derived in essentially the same way. Starting 
with Eq. (p,15|) , imposing Pp = 1 for small momenta, and the general conditions ( |2.1|) , (|2.2| ) 
and (|2.3|), we finally end up with 

Rtk(q) = i (7?-iW 2 -<7 2 ), (2.20) 



normalised with cp = \. In terms of a dimensionless function rp(q 2 /k 2 ), Eq. ( |2.2(J| ) becomes 

r F , opt (v) = (4= - l) ©(1 - V) ■ (2-21) 

and Tf(\) = 1- The non-analyticity of Eq. ( |2.21| ) is a direct consequence of Rp t k having 
only one mass dimension. We shall see below that it is of no harm for the computation 
of fermionic flows because Eq. ( |2.21| ) enters only in specific combinations such that the 
non-analyticity disappears. 



III. DERIVATIVE EXPANSION 

The flow equation is a functional differential equation, which, from a technical 

point of view, is equivalent to infinitely many coupled partial differential equations for the 
couplings parametrising the effective action IV A number of different systematic approxi- 
mation procedures for flows are known. In this section, we consider flows to leading order in 
the derivative expansion, based on expanding the operators of the effective action according 
to the number of derivatives ||29|| . This leads to a closed set of coupled partial differential 



equations for the coefficient functions. We discuss the main structure of optimised flows 
and contrast it with generic ones. 



A. Specific flows 

In order to make our subsequent reasoning more transparent, it is useful to have an 
explicit example at hand. To that end, we consider an 0(A)-symmetric real scalar field 
theory in d dimensions, the linear sigma model. To leading order in the derivative expansion 
we make the Ansatz 

T k = J d d x (u k (p) + \z k {p)d^ a d^ a + ^Y k (p)d^d„p + 0(«9 4 )) (3.1) 

for the effective action, with p = ^0 a a . For Af^l, there are two independent wave func- 
tion factors Z k and Y k beyond leading order in this expansion (cf. Section |V|). To leading 
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order in the derivative expansion, the flow equation ( |1 . 1| ) reduces to a flow for the effective 
potential, dtU k . The main physical applications concern the non-trivial Wilson-Fisher fixed 
point in d — 3 and the computation of related universal quantities. 

Inserting the Ansatz (|3.1| ) into the basic flow equation, and using Z = Y = 1, yields |30 



au _ 1 f A ( (N - l)d t R k (g 2 ) d t R k {q') \ 

1 k 2 J (27TY [q* + R k (qi) + U' k (p) + q* + R k {q>) + U' k (p) + 2pU»{p)) ' 1 ' ] 

It is a second order non-linear partial differential equation. One easily recognizes the contri- 
butions from the iV — 1 "Goldstone" modes and the "radial" mode. A similar flow equation 
has been obtained for the wave function renormalisations Z k and Y k ||30|| . The momentum 
integration is regularised in the UV, owing to the regulator term d t R k {q 2 ) in the numerator, 
and in the IR due to R k (q 2 ) in the denominator. 



B. Generic flows 

For convenience we perform the angular part of the momentum integration and rewrite 
the right-hand side of Eq. (|3.2| ) in terms of so-called threshold functions [|30| as 



= 2MN - DkX + 2 Vd kX {miMM^ . (3.3) 

The constants Vd are given by 

vj 1 = 2 d+1 7r d/2 r(f ) , (3.4) 
and the functions £^(uj) are defined as 

While the flow (|3.3| ) is specific for the theory defined by Eq. (|3~T|) , the functions (|3.5|) are not. 
These functions describe the generic structure of the flow to leading order in the derivative 
expansion. The flows for different indices n > are related by 

dJ d n (co) = -(n + 6 nfi )£ d n+1 (co). (3.6) 

Therefore, it suffices to study the flows £q(uj). 
The fermionic analogue of the flow ( |3.5|) is |24| 

and Eq. (|3.6|) holds equally for Eq. (|3.7|) . Notice the additional factor 1 — r F in the inte- 
grand, which arises due to the Dirac structure of Eq. (|2.20|). We have used the normalisation 
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condition rp{\) = 1. 

It is evident that the characteristics of the flow, determined by the choice of Rk, are en- 
tirely encoded within the functions Eqs. ( |3.5| ) and ( |3.7| ) (or similar functions to higher order 
in the derivative expansion). For a generic regulator, these are complicated functions of the 
fields, which can be computed explicitly only for very specific choices for the regulator. 

Two properties of generic flows given in terms of £q{uj) and £p (uj) are worth being 
mentioned. First of all, from their very definition and the constraints imposed on the 
regulator function, we conclude that any function (|3.5| ) for n = decays at most as 1/u 
for u — > oo |p0| . Therefore, they describe the decoupling of "heavy" modes from the flow, 
which is a manifestation of the decoupling theorem ||19|| . Secondly, all flows have a pole in 
C + uj, where C denotes the gap. Both the analytical structure and the strength of the pole 
depend on the regulator. From the general requirements for regulators, and the explicit 
form of ( |3.5| ), it follows that the pole for n = cannot be stronger than a simple analytical 
pole ~ 1/(C + uS). The pole of threshold functions has important physical implications. 
It determines the approach to a convex effective potential for theories within a phase of 



spontaneous symmetry breaking |31 



C. Optimised flows 

Now we turn our attention to the optimised regulators introduced in Eqs. ( [2.18| ) and 
( p.20| ). The evaluation of Eq. Q3.5|) is particularly simple because the 0-function cuts off 
the momentum integration. Using Eq. ( 2.19| ), Eq. ( |3.5|) reduces to two terms, 



= — t dyyi- 1 + r dyy^ 1 ( y ~ mi - y) . (3i 
1 + wJo yy Jo yy l-yQ(l-y)+w 



In the first term, the momentum integration is cut-off above y < 1. The function P 2 (y) 
remains a constant in this momentum regime, which allowed to move the u;-dependent term 
in front of the momentum integration. The integrand of the second term contains products 
of distributions. Since the integrand is proportional to ~ (1 — y)S(l — y) the second term 
vanishes identically, independently of the specific implementation for the O-function. The 
remaining momentum integration of the first term becomes trivial and gives 

We used the normalisation r(i) = 1 and hence P 2 = 1 for y < 1. 
For fermionic flows (|3.7|) and the regulator (|2.20|) , we find 

£%Ju) = / dyy d/2 ~ l + 2 / dyy2 +1 -XJL — vy — . 

F ' oK } l + oo Jo Jo [^y + (l-^y)&(l-y)} 2 + w 

(3.10) 
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The first term has a restricted momentum integration due to the cut-off provided by the 
G-function. The second term is more involved, and the integrand even contains products 
of distributions. Notice, however, that it contains the factor ~ (-^ — 1)5(1 — y) which 
is proportional to ~ (y — 1)5(1 — y). Therefore, the second term vanishes identically and 
independent of the parametrisation of the distributions and their products. The evaluation 
of the first term gives finally 

and is identical to the bosonic flow. 



D. Discussion 



The flows described by the functions ( ^9|) and ([3.11D have the simplest asymptotic struc- 
ture for uj — > oo. This implies that heavy modes decouple "the fastest" from the flow for 
optimised regulators. For comparison, the sharp cut-off leads only to a logarithmic de- 
coupling ~ \nu. Also, the decoupling does not depend on the particular theory studied 
(i.e. the dimension), unlike the case for polynomial regulators Rk ~ q 2 (k 2 /q 2 ) b . Further- 
more, the flow described by the functions ( |3.9| ) and fl3.11| ) have the simplest and strongest 
pole structure for C + uj — > + . The pole is a simple analytic one, which is not the case 
for generic regulator functions. An immediate implication of this structure is that the op- 
timised flows (|3.9|) and fl3.11|) lead to a logarithmically smooth approach towards a convex 



effective potential. This is very different from the sharp cut-off case, where the approach is 
only exponentially ||. A detailed presentation of these results is given elsewhere. 

For completeness we quote the results for the flows ( |3.9| ) and ( |3.11| ) for arbitrary normal- 
isation. While the normalisation is of no relevance for a theory containing only bosonic or 
fermionic degrees of freedom, their relative normalisation can become important for theories 
containing bosons and fermions. The normalisation conditions ( |2.5|) and ( 2.14[) correspond 



to r(ce) = 1 and vf{cf) = 1, which can always be imposed because the functions r(y) and 
rp{y) range between < r,r F < oo. The optimised gaps are C opt = 2c b and Cp opt = 4cp. 
For arbitrary eg the flow Eq. (|3.9| ) is obtained as 

= (3-12) 



In the fermionic case we find 



d 2cb + uj 



* oM = (3,3) 



d Acf + uj 



for the rescaled analogue of Eq. (|3.11|) . 
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Finally, we note that flows £q(uj) ~ have been used earlier in the literature p2^p4 
however, without the explicit knowledge of the corresponding regulator function. These 
trial functions are sufficiently simple to allow for analytical considerations. The motivation 
for their use was based on the observation that the generic threshold function Eq. (|3.5| ) 
decays at most as w" 1 for large u. This suggested that a regulator may exist which leads to 
£q(w) = Ad(C + cj) -1 . Let us show how the normalisation can be fixed from consistency 
arguments. We use the universal relation ^ n (0) = 1 P0|| , which holds for d = 2n dimensions, 
to identify the prefactor as A 2n = ^C n+1 . The analytic continuation to arbitrary dimen- 
sions leads finally to our results (|3.12|) and (|3.13| ). This reasoning shows that the Ansatz 
£q(co>) = 2C d / 2+1 (C + uj)~ 1 is self-consistent. However, we rush to add that these consistency 
arguments are necessary conditions, but not sufficient ones: only the explicit form of the 
regulator - as given by Eqs. ( [2.18| ) and ( |2.20| ) - finally justifies the few earlier computations. 
In addition, Eqs. (|2.18|) and (|2.20|) are explicitly required for the computation of the flow 
at finite temperature (see Section |TV| ) or to higher order (see Section [V]). 



IV. THERMAL FLUCTUATIONS 



In this section we apply our reasoning in the context of a quantum field theory coupled 
to a heat bath at temperature T, and to leading order in the derivative expansion. We show 
that optimised flows, as opposed to generic ones, disentangle the different contributions 
related to thermal and quantum fluctuations, respectively. These properties are realised, on 
the level of the flow equation, in terms of an important factorisation. This leads to better 
convergence properties of the flow itself. Approximate solutions of the flow correspond to 
better approximations of the physical theory. 



A. Imaginary time formalism 

To be explicit, we consider a bosonic or fermionic field theory at thermal equilibrium at 
the temperature T within the Matsubara formalism. This implies that periodic (resp. an- 
tiperiodic) boundary conditions for the bosonic (resp. fermionic) fields are employed. As a 
consequence, the go integration in the flow equation (|1.1| ) is replaced by a sum over Mat- 
subara modes m = 0, ±1, ±2, . . .. The trace in Eq. (|1.1|) contains a momentum integration, 
which is then substituted as 



/ 



^ TZfjt^L. (4,) 



(2n) d (2ttj 
In the integrand of Eq. (|1 . 1| ) the go variable is replaced by 

go -> 2nc m T (4.2) 

where 
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c m = m for bosons (4.3) 

c m = m + | for fermions . (4-4) 

It is also useful to introduce the variable 



T 



2nT/k (4.5) 



for the following considerations. The replacement ( (4.2| ) implies that functions turn 
into temperature dependent functions £^(u,t). We show that this function factorises for 
the regulators (|218|) and ( ^20|) . 



B. Dimensional reduction and fermion decoupling 

Let us review a few basic facts known for generic flows at finite temperature within the 



imaginary time formalism |3(| , |35| , |24| . |36 



Bosonic fields within the Matsubara formalism display the phenomenon of dimensional 
reduction at high temperature. This means that for T large enough all non-vanishing 
Matsubara modes are suppressed due to effective masses ~ mT for the Matsubara modes 
with m 7^ 0. Only the m = mode survives in this limit, leading to an effective theory in 
(d — l)-dimensions. For a generic bosonic regulator, the finite temperature flow is given as 

v d 2ir ^ Jo P 2 (y + c m r 2 ) + u 

The function P 2 is defined in Eq. (|2.9|) . The asymptotic regime where only the m = 
Matsubara mode contributes is reached for r — >• 00. From Eqs. ( |4.6| ) and ( |3.6| ), we deduce 

On the other hand, the limit r — >• eventually switches on all higher order Matsubara 
modes. It is straightforward to verify that 

e d n (u,r^0) = £ d n (u). (4.8) 

The asymptotic limits for T — > 00, Eq. (|4.7|), and T — > 0, Eq. ([4.8|), display dimensional 
reduction for bosons as a function of temperature for generic regulator function |[35|| . 

Fermions at finite temperature within the Masubara formalism can be treated in es- 
sentially the same way. However, they happen to have no m = mode as antiperiodic 
boundary conditions have to be used on the go-integration. Hence, fermions do not display 
dimension reduction. Rather, they decouple completely from the RG flow once the small- 
est Matsubara mode is larger than the scale k. These properties can be read-off from the 
temperature-dependent flow. For a generic fermionic regulator, the flow £p (uj) at finite 
temperature is defined as 
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The function Pp is given in Eq. (|2.15 ). The asymptotic regime where the fermions decouple 
completely is reached for r — > oo. From Eq. (|4.9|), we deduce that 



4„(u;,t->oo) = 0. (4.10) 
Again, the limit r — > eventually switches on all higher order Matsubara modes such that 

4>,r^o) = 4». (4.ii) 

The asymptotic limits ( |4.10| ) and ( [4.1 1| ) describe the decoupling of fermions in the high 
temperature limit for arbitrary dimension and generic regulator function. 



C. Optimised thermal flows and factorisation 

We now turn to the optimised regulators ( [2.1S| ) and ( |2.20| ). For this case, the flow 
can be computed explicitly. Inserting Eq. (|2.18|) into Eq. ( |4.6|) , and following a reasoning 
analogous to the one after Eq. (|3.8|), we find 



ti{u,T) = B d {r)l d n {u) (4.12) 
with the temperature dependent function 

Notice that the temperature effects have factorised. This implies that temperature cuts off 
all amplitudes u in the same manner. This is not the case for a generic regulator. 

Let us discuss the thermal threshold factor B^t). In Fig. 2 the thermal threshold factor 
Bd(r) is displayed for d — 2,3 and 4 dimensions. Every single Matsubara mode contributes 
to Eq. ([4.13|) proportional to 

-rfl-cM^efl-c^r 2 ) . (4.14) 



The B-function is a remnant of the regulator ( |2.19| ) and cuts the m-th Matsubara mode 



off as soon as k < c m T/27c. The factor r stems from the go-integration and the factor 
(1 — c^r 2 )^" 1 )/ 2 from the d — 1 dimensional integration over spatial loop momenta |q|. 
These functions vanish outside the interval < r < l/c m . At the upper end they behave 
like (l/c m — t)^ 1 ^ 2 and vanish linearly with r at the lower end. This structure explains 
the spikes observed in Fig. 2, which are located precisely at the points r = l/c m and due 
to the decoupling of the ±c m -th Matsubara modes. Indeed, for r > 1 only the m = 



Matsubara mode yields a contribution to B^t) in Eq. (|4.13j). The asymptotic regime 
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where only the m = Matsubara mode contributes is reached already for r > 1 with 
Bd{r > 1) = Tdv d -i/27rv d (d - 1), or 



v d -i T 
v d k 



id-l, 



(4.15) 



Notice the difference to Eq. (|4~7| )- Decreasing r below r = l/c m eventually switches on 
the ±c m Matsubara modes. For r close to the points l/c m , the term ( [4.14| ) increases 
as (l/c m — r)( d_1 )/ 2 for decreasing r. This power law explains why the spikes are more 
pronounced in lower dimensions.^] In the limit r — > it is straightforward to verify that 
Bd(r — > 0) — > 1 which implies 



(4.16) 



This asymptotic limit is the same as Eq. 




0.6\ 



0.5 1 1.5 

T = lixTjk 



Figure 2: Dimensional reduction for bosons, described by the bosonic thermal function 
B d (T) defined in Eq. ( 4.13| ). d = 4: full line, d = 3: dashed line, d = 2: dashed-dotted 
line. 



**To higher order in the derivative expansion, the spikes are smoothed-out for non-trivial wave 



function renormalisation, cf. Sect.VC 
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Turning to the optimal fermionic regulator ( |2.21| ), the momentum integration in Eq. ( |4.9| ) 
can be performed explicitly to give 



i d F (u,T)=F d (T)i d F (u) . 



(4.17) 



As in the bosonic case, the temperature effects factorise from the threshold effects. The 
fermionic thermal threshold factor F<i(t) is given by 



FJt) 



d v d _i r 



2 2 
C T 



d-1 
2 



e(i- 



2 2 



(4.18) 



d — 1 Vd 2n 

Eq. ( 4 . 1 8|) is identical to its bosonic counterpart ( |4.13| ) except for the Matsubara sum which 
runs over c TO = ±|, ±|, . . . in Eq. ft4.18|) . In Fig. 3 we have displayed the function F d {r) for 
d = 2,3 and 4 dimensions. Again, the spikes have the same origin as in the bosonic case 
and the same reasoning applies. The high temperature limit at which the fermions decouple 
completely, is already reached for k < nT, 

t% n {u, r > 2) = . (4.19) 

Notice the important difference to Eq. ( f4.1U| ), where the decoupling of fermions is only 



asymptotic. The limit r — > is equivalent to Eq. (4.11 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 








Fd(r) 




• » \ * 
• \ \ 

\ * ' 









0.5 



2.5 



3 



1 1.5 2 

t = lixTjk 

Figure 3: Fermion decoupling at finite temperature described by the fermionic thermal 
function F d (r) given by Eq. ( 4.1 8| ) . d = 4: full line, d = 3: dashed line, d = 2: dashed- 
dotted line. 
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D. Discussion 



The optimised regulators (|2.18|) and Q2.20|) correctly describe dimensional reduction and 
fermion decoupling. In addition, they lead to a thermal factorisation of the flow as observed 
in Eqs. (|4.12j ) and ( 4.171) . From a physical point of view, this fact is easily understood. 
The imaginary time formalism compactifies the time direction and the temperature mod- 
ifies the temporal momentum modes of the fields. The corresponding Matsubara mode, 
when compared to the infrared scale k, leads to a thermal decoupling. To leading order 
in the derivative expansion, the optimised regulator makes the temperature blind for the 
quantum fluctuations. It cannot distinguish between amplitudes of constant fields. In turn, 
the quantum fluctuations are sensitive to the field amplitudes, which are responsible for the 
mass decoupling, similar to the case for vanishing temperature. Clearly, these two effects 
have different physical origins. As a consequence, it is natural to employ a regulator which 
reflects this factorisation on the level of the flow equation. 



For a generic regulator function, the flows (^6) and (|4.9| ) are complicated functions of 
both the temperature and the field amplitudes. They reflect dimensional reduction and 
fermion decoupling. Typically, however, they do not factorise. This simply means that a 
generic ERG flow entangles thermal and quantum fluctuations even to leading order in the 
derivative expansion. This is a direct consequence of the regulator term, whose coupling 
to the different operators in the effective action is leading to a field-dependent thermal de- 
coupling of the different modes on the level of the flow equation. This entanglement is of 
no relevance if the flow can be solved exactly. In turn, for an approximate solution of the 
flow, the factorisation on the level of the flow equation is most helpful. It avoids a mixing 
of thermal and quantum fluctuations in a regime where they can be disentangled, thereby 
minimising possible artifacts due to the specific regulator function. As a consequence, the 
flow itself is stabilised, and expansions of the flow show much better convergence behaviour 
towards the physical theory. More generally, it is expected that this line of reasoning applies 
for generic optimised regulators. 

Finally, the factorisation is very helpful for numerical solutions of flow equations. In 
the generic case, one two-parameter function has to be fitted in order to describe the flows 
( |4.6p or ( }4.9| ). In turn, only two one-parameter functions are needed once they factorise as 
in Eqs. ( |4.12|) and ( f4.17| ). This simplification is substantial, and even more so, because the 
functions (}4.12|) and ( |4.17| ) have a very simple analytical form. 



V. DERIVATIVE EXPANSION TO SECOND ORDER 

In this section, we apply the optimised regulator to higher order in the derivative expan- 
sion. We first discuss the general structure of the equations. Furthermore, we show that a 
simple factorisation of the flow takes place for field-independent wave function renormali- 
sation. The physical origin of the factorisation is discussed, and its realisation on the level 
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of the flow equation leads, as in the thermal case, to better convergence properties of the 
flow and the derivative expansion. For technical details on the computations, we defer to 
the Appendices . 



A. Wave-function renormalisation 



In the preceding sections we have restricted the discussion to the leading order in a 
derivative expansion. This implied the vanishing of the anomalous dimensions rj^ = rj^ = 
or Ztj, = = 1. To higher order in the derivative expansion the multiplicative renormali- 
sation of the fields has to be taken into account, according to 



Once higher derivative terms are included in the Ansatz for the effective action, additional 
flow equations for the corresponding coefficient functions like Z^kiq 2 , p) and Z^^(q 2 , p) have 
to be studied. The wave function renormalisations are functions of the scale parameter k 
and can depend as well on momenta q 2 or on the mean fields p. To second order in the 
derivative expansion the wave function renormalisation is evaluated at a particular momen- 
tum scale q 2 = k 2 which fixes the renormalisation conditions. Typical choices for ko are 
k = and k = k. 

The most important new ingredient at this order is the scale- and field- dependence of the 
wave function renormalisations. In the example defined through Eq. (|3.1|) , these are given 
by the functions Zk and Y^. Here, the function Zk is responsible for the renormalisation of 
the N — 1 "Goldstone" modes, which differs from Zk = Zk + pYk for the "radial" mode. The 
fact that different wave function renormalisations appear to second order in the derivative 
expansion (depending of the theory considered) is of no relevance for the following discus- 
sion of the flows. The parametric dependence of the flow on either Zk or Zk is the same. 

Let us introduce an additional function Zk(p) as 



We have factored out a constant term Zk(po) chosen at an arbitrary reference point. We 
have chosen the reference point p = po which fixes the renormalisation of the fields for all 
momenta. Typical choices are either po = 0, or po — the minimum of the scale-dependent 
potential. The split ( |5.3|) allows to separate the non-trivial field-dependence, contained in 
Zk(p), from an overall renormalisation contained in Zk{po)- The factor Zk is normalised as 

Zk(po) = 1. 

In order to provide a simple form for the flow it is useful to introduce the field- 
independent factor Zk(po) into the regulator function, 



4> - O 




(5.2) 




(5.3) 
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Rk — > Z^^Rk 
RF,k Zf,kRF,k 



(5.4) 
(5.5) 



The flow equation, when written in terms of renormalised variables ( |5.1| ) and ( |5.2j ), receives 
additional contributions proportional to the anomalous dimensions 

= ~9t In Z^ k , (5.6) 
Vi> = ~9t In Z^ k , (5.7) 

because the derivative d t Rk in the flow equation now acts also on the explicit scale- 
dependence contained in Z k . 



B. Generic flows to second order 

To second order in the derivative expansion the flow has turned into a function of the 
field amplitudes ou, the anomalous dimension rj and the field dependent function z. The 
corresponding bosonic flow is defined as 



£%{u;,z, V ) = dyy2- 1 * ^ 2 / * . (5.8) 
Jo y[z + r(y)\ + u 

Notice that the pole structure of the flow is changed, owing to the function z(p). The 
effective inverse propagator becomes a function of the fields: 

P 2 (p,y)=y[z(p)+r(y)]. (5.9) 



The location of the pole of Eq. ( p.8| ) at — lu = C(p) = min^o P (p, y) has turned into a func- 
tion of the fields. For the optimised regulator (|2.19|) the pole is located at C = min{l, z}. 
Compared to the leading order in the derivative expansion, the pole structure is modified 
once z ^ 1. 

For the fermionic case, the flows are given as 

Jo y\z F + r F (y)\ 2 + lu 

The pole structure changed as well, as follows from 

P 2 (p,y)=y[z F (p)+r F (y)} 2 . (5.11) 

The location of the pole of Eq. (|5.8| ) at — uj = C F (p) = min y > P F (p, y) has turned into a 



function of the fields. For the regulator (2.21) the pole is located at C F = min{l, z F }. The 



pole structure is modified compared to the leading order in the derivative expansion once 
z F ^l. 
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C. Optimised flows to second order 



We now turn to the optimised flows and discuss their structure at second order in the 
derivative expansion. We refer to the Appendices for all technical details. 

In the case of a generic wave-function renormalisation where Zk(p) is a non-trivial func- 
tion of the fields, the function £q(lu, z,rj) as defined in Eq. ( |5.8| ) can be evaluated explicitly 
for the regulator ( [2.18j ). The structure of the flow is as follows. Consider the denominator of 
Eq. (|5.8| ), given by y(z + r) +u. It can be rewritten as 1 + u + y(z — 1) + (y(l + r) — 1). The 
last term in brackets vanishes for the optimised regulator fl2.18|) because the integration is 
restricted to y < 1. The remaining term can be written as the product (1 + u)(l — jr^y)- 
Notice also that the numerator of Eq. (|5.8|) neither depends on z nor on u. These obser- 
vations lead to the following conclusions. First, and apart from an overall cu-dependence 
~ (1 + uj)^ 1 , the optimised flow depends on z and u only through the variable 

1 — z , 
£ = — — • 5.12 

1 + UJ 

Second, the optimised flow factorises into a leading order term (|3.9|) and a remaining factor 

£ d (cu,^ V ) = £ d (cu)B d ^,r ] ) , (5.13) 



because the denominator of Eq. (|5.8|) contains a momentum-independent factor (1+lu). Here 
we have introduced £ d (u, £(a>, z), rj) = £q(uj,z,t]). It it interesting to see that the structure 
of the optimised flow is still quite simple. An integral representation of Ba(C,,i]) is given in 
Eq. ( |B.3| ). (For all £ < 1, Bd(£,r)) can be expressed in terms of hypergeometric functions, 
cf. Eqs. (P-4j ), (p.5| ) and ( |B.6|) ; closed expressions of Eq. ( |5.13|) for d = 4, 3 and 2 dimensions 
are given in Eqs. ( |B.7| ), ( |B.8| ) and ( |B.9| ), respectively.) For |£| < 1 the function Bd(£,rj) can 



be Taylor-expanded in arbitrary dimensions, to wit 



Bi(Ul) = SrTf I 1 " J^Tn) = l ~J+2 + tti +oie ^ ■ (5 ' 14) 



The series representation (^.14 ) is best suited for the flow as long as |£| remains small. 



This corresponds to either the limit of a field-independent wave function renormalisation 
z(p) = 1, or, for any z, to the limit of large amplitudes u. From Eq. Q5.14 ), we obtain for 
Eq. ( [5.1 3| ) to zeroth order in £ 



«^ )= «. Hl -^j +0 (^,jy . (5.i5) 

We note that the ?7-dependent correction in Eq. ( |5.15|) has the same functional dependence 
on the amplitude u as Eq. (p.9|). Stated differently, the optimised regulator leads to a sim- 
ple factorisation in both the decoupling limit uj ^> 1 and for the case of a field-independent 
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anomalous dimension £ = 0. 

For completeness we give also the result for the bosonic flow at finite temperature. The 
corresponding flow £q(uj, £, r, 77) still factorises as 

%(u,£,T,Tl) = B d &T,Tl)el(u). (5.16) 

This is the generalisation of Eq. (|4.12|) to second order in the derivative expansion. We only 
have to replace the function 77) by its temperature-dependent counterpart B d (£, r, 77) 
[cf. Eq. ( |A.5|) ]. It is straightforward, if tedious, to establish explicitly that the function 



B d (^,r,r]) represents dimensional reduction in precisely the same way as B d {r) = B d (^ = 
0,r, 7] = 0). Let us consider the most interesting case, which is the leading order in £ <C 1. 
In this limit, Eq. ( |5.16j ) reads 

B d ^,r,r ] ) = B d (r,r ] ) + 0^^r ] ) . (5.17) 

The function B d (r, r/) can be expressed as 



BAr,v) = B d (r)- V ^ V -^I-Y,^-c 2 m T>) 2 e(l-4r 2 ). (5.18) 

a 1 U d Z7T m 

The function B d (r), Eq. (|4.13|) , has been discussed in Section [IV]. The new ingredient, 
beyond leading order, are given by the corrections ~ 77 in Eq. ( |5.18| ). Every single Matsubara 
mode contributes as 

-rfl-c^^efl-c^ 2 ) . (5.19) 



d+l 



Compared to the leading order contributions ( f4.14| ), we notice that Eq. ( |5.19| ) follows from 



Eq. (|4.14j) for d — > d + 2. The reason is very simple. In the flow equation, the anomalous 
dimension is proportional to a term containing an additional factor ~ q 2 , which effectively 
increases the momentum measure by two dimensions. This has an immediate consequence. 
The thermal decoupling in Eq. ( |5.18| ) proportional to ~ 1] is much smoother than the leading 
order decoupling, simply because the spikes are less pronounced the higher the dimension. 
Therefore, the spikes observed in Fig. 2 are smoothed-out once 77 (and £) are non- vanishing. 

In the opposite regime where j^j- 3> 1, only a few leading terms of the series ( |5.14| ) 
have to be retained. This limit is of relevance close to the pole region of the flows lo — > — 1, 
or in the region of large z 3> 1. From the explicitly resummed expressions ( p.7|) , (|B.8|) and 
( |B.9| ), we conclude that a factorisation as 

e*(u,z,ri) = f d (u,z)(l-^ (5.20) 

holds true, and f 4 (z) = (z - f 3 (z) = 2f 4 (z) and f 2 (u, z) = f 4 (z) ln(f±^). 

It is not surprising that a similar structure is found for fermionic flows. The correction 
term due to the substitution Eq. ( p.4| ) simplifies Eq. ( |5.10 ) to 
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*M ,„ „) = f iyy^ 11 + f; F ~ 1)111 T;;:' 1 -^" . (5.21) 
Jo [I + y/y{z F - + w 

Eq. (|5.21 ) factorises as 

£ d Ffi (uJ, z F , rj^) = £ d Ffi (uj) F d (u, z F , rty) . (5.22) 

The function F^lu, z,rj) can be expressed in terms of hypergeometric functions. At finite 
temperature, and for z = 1 and r] = 0, it reduces to Eq. Q4.18Q . Here, we are only interested 
in the structure of the flow for a nearly field-independent wave function renormalisation, 
z « 1, or for the decoupling limit. We find 

FM ,„ „, = 1 _ JM_ + (l _ _*_) (l _ + . . . ( , 23) 

The two leading terms in Eq. ( |5.23|) show that also fermionic flows factorise for field- 
independent wave function renormalisation. 



D. Discussion 



The structure of the flow has increased to second order in the derivative expansion. Let 
us discuss first the case of a field- independent wave function renormalisation z = 1. The 
corresponding flows ( |5.15| ) and (|5.23|) for the optimised regulators factorise, similar to the 
thermal case to leading order in the derivative expansion. Physically speaking, this structure 
can be made plausible as follows. The flow, when written in terms of the renormalised fields 

- and under the assumption that the renormalisation is momentum- and field-independent 

- depends, in addition to the fields, only on the anomalous dimension. The anomalous 
dimension is field independent, and, in consequence, unable to distinguish between fields of 
different amplitudes contained in u, which parametrise the quantum fluctuations. There- 
fore, it is natural that the flow factorises the contributions induced through 77 from those 
induced by the amplitudes u>. The disentanglement is realised by the optimised regulators.^] 

In turn, a generic flow does not reflect this factorisation. Rather, it leads to an entangle- 
ment between the renormalisation of the effective potential induced by the infrared regulator, 
and the renormalisation parametrised by a field-independent anomalous dimension. This 
is immediately evident from the observation that the ^-dependent and the ^-independent 
contributions to the flow of the effective potential have different functional forms as func- 
tions of the fields. At this level, the entanglement is due to the regulator, which modifies the 



ttFrom the definition of Eqs. (]^J^) and (|5.10| ) it follows that all homogeneous regulators with 
r (y) ~ y r '(y) ( or r F{y) ~ y r ' F (y)i respectively) factorise the anomalous dimension from the field- 
dependent part of the flow. 
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coupling amongst all operators of the effective action. As mentioned in the thermal case, the 
entanglement is of no importance for the full solution to the flow. In turn, the factorisation 
is very useful for approximate solutions. It leads to more stable flows because irrelevant 
couplings, entirely due to the regulator, are removed. The same reasoning as given at the 
end of the previous section applies. 

For the thermal bosonic flow ( {5.16Q , we notice that the dependence on the anomalous 
dimension enters the thermal factor Bd(T,r]). In particular, the thermal corrections do not 
factorise from those due to a field-independent anomalous dimension. This structure can 
be understood as follows. The wave function renormalisation enters the momentum trace 
as a multiplicative renormalisation proportional to the kinetic term q 2 . At finite tempera- 
ture within the imaginary time formalism, the spatial and the temporal loop momenta are 
treated in an unequal way. Hence, thermal fluctuations couple in a non-trivial manner to 
the anomalous dimension of the fields. This implies that the temperature-dependent factor 
itself is modified due to the anomalous dimension, which provides the physical reason way 
no factorisation of the temperature effects from the anomalous dimension are expected in 
the first place. 

For the case of a field-dependent wave-function renormalisation, a simple factorisation 
similar to Eq. (|5.15| ) is not expected, simply because the wave function renormalisation is 
a function of the fields. Hence, the wave function renormalisation can distinguish different 
field amplitudes, in contrast to the field-independent case. However, two observations are 
still worth mentioning. First of all, we observe a partial factorisation, which is evident from 
Eqs. ( |5.13| ) and ( |5.16| ). This structure is based on the fact that the ^-dependence enters 
only through the variable (|5.12|) , as opposed to the generic case. For •fc^l <C 1, only a few 
leading terms have to be retained from the explicit series ( |5.14|) . It follows that each power 
of Jjqj^ is renormalised proportional to the anomalous dimension and an order-dependent 
numerical coefficient. Secondly, the limit for fc^ ^> 1 again allows for a simple factori- 
sation, as follows from Eq. (|5.20|) . Here, the wave function renormalisation can no longer 
distinguish field amplitudes, allowing for this simple structure. 

A final comment concerns the numerical prefactors ~ rj as found in Eqs. ( |5.15| ) and ( |5.23| ). 
We emphasise that the coupling of the anomalous dimensions to the effective potential is, 
apart from the field dependence, dimensionally suppressed - by factors for bosons and 

for fermions - as opposed to the leading order contributions. This additional suppression 
is noteworthy because the convergence of the derivative expansion is controlled by small 
anomalous dimensions of the fields. Here, we have just shown that an expansion performed 
with an optimised regulator leads to an additional dimensional suppression of the back- 
coupling of the anomalous dimension to the effective potential. A more detailed discussion 
of this observation will be given elsewhere. 
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VI. PROPER-TIME REGULARISATION 



In this section we leave aside the conceptual framework of the ERG based on a 
momentum-scale regularisation and address flows based on an operator cut-off regulari- 
sation. Our aim is to provide the analog of the optimised regulator ( 2.18 ) within the 



proper-time regularisation method. For a more detailed comparison with the exact renor- 
malisation group, we refer the reader to Ref. |2l| . 

A simple flow has been derived from a one-loop expression for the effective action which 
is UV and IR regularised using a Schwinger proper-time representation of the operator trace 



f3T|j , amended by a regulator function (A, s) within the proper-time integral [|S8]| . The 
flow with respect to the infrared scale parameter k follows from a 1-loop improvement as 



ti k 



(A, 



Tr exp (sT^ 



(6.1) 



We refer to this flow as the "proper-time renormalisation group" (PTRG). It describes 
the partial resummation of perturbative diagrams. The proper-time regulator function 
plays the role of the momentum regulator within the ERG. The flow ( |6.1|) is governed 
by the IR scale k. Following Ref. pD[ , we introduce a dimensionless function f(x) as 
fj, (A, s) = f(A 2 s) — f(k 2 s) and require f(x — > oo) = 1 and f(x —> 0) = 0. This ensures 
that the usual Schwinger proper time representation is reached in the UV limit. 



We are not aware of a simple and generic optimisation criterion, analogous to Eq. (|2.11| ), 
which derives from within the PTRG formalism. Furthermore, the flow (|6.1|) has no path 



integral derivation, which makes a conceptual reasoning much more difficult. However, it 
is still possible to show that a function f op t(x) exists which is equivalent to the optimised 
ERG regulator (|2.18|) to the leading order in the derivative expansion. 

To that end, we apply Eq. (|6.1| ) to a iV-component real scalar theory in d dimensions 
and to leading order in the derivative expansion. Using the Ansatz (|3.1|) the flow for the 
effective potential Uk(p) with p = ^4> a 4>a becomes 



d t U k (p) = !(47r)-3 j o -^a/«>(A, a 



(u' k (p)+2pU»(p)) 



+ {N-l)e 



sU'(p) 



(6.2) 



This flow is identical in form to the ERG flow ( |3.3| ), if we replace the ERG flow in Eq. ( |3.3| ) 
by the proper-time flow 



1 — 

dxx 2 (d t f{x)) exp(— xuj) . 



(6.3) 



Here, the integration variable is x = k 2 s and stems from the proper-time integration, in 
contrast to Eq. ( p.5|) , where y = q 2 /k 2 stems from the momentum trace. Now, consider a 
specific class of proper-time regulator functions: 
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m = , dj(x) = (6.4) 

1 [m) 1 (m) 

We have introduced a free parameter m describing different regulators, and the incomplete 
T-function T(m,x) = Jq dyy m ~ x e~ y . This yields the simple expression 

^M = r(m r 7f (i) (l + ^- m (6-5) 
f (m) 

which agrees with Eq. (|3.9|) for m — 1 + |, or 

The optimised proper-time regulator (|6.6| ) corresponds to the optimised regulator (|2.18|) 
within the ERG approach. Hence, it is possible to identify an optimal regulator function 
for proper-time flows, owing to their close similarity to the ERG to leading order in the 
derivative expansion. 



VII. CONCLUSIONS AND OUTLOOK 

This study was motivated by two observations. First, an application of the ERG to 
realistic physical problems is bound to certain approximations. Second, approximate so- 
lutions of flow equations depend spuriously on the infrared regulator. Combining these 
observations, it became obvious that an understanding of the spurious scheme dependence 
is mandatory in order to provide predictive power for approximate solutions. Previously, 
we showed that the gap of the full inverse propagator controls convergence properties of 
approximate solutions ||15|| . It has also been shown, based on the computation of critical 



exponents for the Ising universality class, that the convergence of the derivative expansion 
is controlled by the gap |L8| . These observations lead to the conclusion that the freedom in 
the choice for the IR regulator can be used to maximise the physical information contained 
within a given approximation or truncation. 

An interpretation of the interplay between the RS function on one side, and conver- 
gence of approximate flows on the other, is as follows. The IR regulator - by regulating the 
flow - modifies the interactions at intermediate scales fc ^ amongst all operators of the 
theory. Eventually, these cancel out for the integrated full flow, but not for approximated 
ones. Hence, changing the RS function for an approximated flow modifies some remaining 
RS dependent terms which cannot be cancelled due to the missing contributions from ne- 
glected operators. Therefore, a "fine-tuning" of the RS function allows to partly incorporate 
higher-order effects within the lower orders of a given approximation. This corresponds to 
an optimisation. 

The present derivation of optimised ERG flows had two ingredients. First, we made use 
of a generic optimisation criterion for bosonic and fermionic fields ||15|| , which states that the 
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gap of the full inverse propagator, as a function of momenta, should be as large as possible. 
This way, the ERG flow is the least singular, and approximations to such flows are expected 
to be much more stable, leading to improved convergence of expansions. Second, we added 
the specific requirement that the effective inverse propagator be momentum-independent in 
the IR region (cf. Fig. 1). 

We have studied specific optimised ERG flows for bosonic or fermionic theories up to 
second order in the derivative expansion and at vanishing and non-vanishing temperature. 
Their specific properties have been discussed in detail at the end of the corresponding sec- 
tions. Here, let us only mention the perhaps most surprising property of optimised flows, 
which is the disentanglement of thermal and quantum fluctuations to leading order in the 
derivative expansion. A similar factorisation occurs for field-independent wave function 
renormalisations. 

More generally, optimised ERG flows owe their main properties to the "flatness" of the 
effective inverse propagator, which extends over the entire momentum region q 2 < k 2 for 
the specific regulator studied here. Other regulators can lead to similar factorisation and 
convergence properties. Prime candidates are given by solutions to the optimisation crite- 
rion: as is evident from Fig. 1, they automatically lead to flat effective propagators - at 
least within a small region about the minimum of the effective inverse propagator. If this 
region extends over the domain where the flow equation receives its main contributions, we 
expect to find equally good flows. 

An important conclusion is that the optimisation ideas discussed here should be useful 
for high-precision computations based on this formalism. Increasing the precision normally 
implies a full computation at the following order of the expansion. Here, we argued that the 
physical results can be improved already within a fixed order of the expansion. Immediate 
applications of optimised flows concern the computation of universal critical exponents for 
TV-component scalar theories in three dimensions, or the study of convex effective poten- 
tials within a phase with spontaneous symmetry breaking. On the conceptual side, it is 
possible to show that the optimisation criterion can be interpreted as a natural minimum 
sensitivity condition, somewhat similar to the principle of minimum sensitivity as employed 
within perturbative QCD. We will leave a detailed discussion of these results to a future 
publication |18[ . 



Our analysis can be extended in a number of directions. For gauge theories, modified 
Ward or BRST identities ensure the gauge invariance of physical Green functions p9| , and 



the optimisation criterion is compatible with such additional constraints. This optimisa- 
tion can also be implemented for field theories at finite temperature within the real-time 
formalism ||36fl . While our present analysis is based on the derivative expansion, it seems 



worthwhile to study optimisations for other systematic expansions like expansions in powers 
of the fields. Finally, it would be interesting to see how these ideas apply to flows at finite 
density or to Hamiltonian flows 
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Appendix A. FLOWS TO SECOND ORDER IN THE DERIVATIVE EXPANSION 

In this Appendix, we derive explicit expressions for the optimised flow for the effective 
potentials to second order in the derivative expansion at both vanishing and non-vanishing 
temperature. 

Our starting point is the flow for the effective potential to second order in the derivative 
expansion. We consider the bosonic flow at finite temperature. All information about the 
flow is parametrised by 

e^, z , T , n) = ^ti— £ t^ dy ^-h^-y-f*W-^ (A.i) 

IA " v d 2-k „ 7o (z-1)(j/ + <£t 2 ) + 1+u> v ' 

Here, u;(p) is the field variable, z(p) a field dependent wave function renormalisation, 77 the 
anomalous dimension, r = 2-nTjk the rescaled dimensionless temperature, and c n = n the 
Matsubara modes in the bosonic case, Eq. ([4.3|). The constants Vd are defined in Eq. (|3.4j) . 

The leading order behaviour of Eq. (|A.1|) is given by the function 

^) = ~. (A.2) 

which follows from Eq. ( |A.1| ) for z — 1, r = and 77 = 0. Factorizing the main building 
block Eq. (|A.2|) from Eq. ( |A.1|) , we notice that the remaining factor depends on both u and 
the variable z only through the combination 

(A.3) 

Therefore, it is most natural to make the variable transform Eq. ( |A3| ) by writing 
£q(uj, z, t, T)) = £q(uj, u>), t, T]), and to rewrite the flow ( |A.1| ) as 



e%(u,Z,T,7i)=$(u)B i (£,T,Ti) , (A.4) 

where 

m,T,v) = ^-X r^dyyl-i 1 - ' <*,- y ~ t 2) 6(1 - 4r 2 ) (A.5) 

dV ^' ' /; 2 v d 2n Z ^Jo yy l-£(y + c 2 n T 2 ) v n ! K ' 

Below, if not stated otherwise, we adopt a simplified notation: functions are evaluated at 
the points £ = 0, z — 1, r = or 77 = 0, if the corresponding arguments are not displayed. 
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With these definitions at hand, we can face the explicit computation of Eq. ( |A.1|) . 

Let us compute the function B d (C,,T,r]) more explicitly. Since the anomalous dimension 
enters only linear in Eq. ( |A.1| ), it is helpful to rewrite Eq. ( |A.5| ) as 



B d {Z,T,r)) = B d (£,T)-riB d (£,r) 



(A.6) 



For £ < 1, the remaining integration over the momentum variable in Eq. ([A. 5 ) can be 
performed. This leads to 



1 7 E^^)(l-4r 2 )^©(l-4r 2 )x 



G? — 1 f d 27T 



The function 



2 F 1 (a,b; c; z) 



r(c) 



r(6)r(c-6) 



[ dt t b - 1 (i-ty- b - 1 (i-tz)- a 

Jo 



(A.7) 



(Ai 



with 2 Fi (a, 6; c; 2) = 2-^1 (&> a; c; 2) denotes the hypergeometric function We also intro- 
duced the thermal amplitude function 



A„(£,r) 



2^2 



1 - <T 



2^2 ■ 



(A.9) 



The index n corresponds to the Matsubara mode. The factors A n only appear in combination 
with the factor 0(1 — ctr 2 \. At the limits, 



A n (£,0)=A(£,r) = l 
A n (e, c ; 1 ) = o. 



(A.IO) 
(A.ll) 



In the same way, we find from Eq. ( A. 5 ) for the term ~ 77 in Eq. (|A.4j) the explicit expression 



2jPi (1, <*±I ; ^ n (£, r)) - g a F a (l, ^; ^; ^4 n (£, 



(A.12) 



Combining Eq. ( |A.7| ) with Eq. (A.12) gives Eq. ( |A.6| ) and hence Eq. ( |A.4j ) explicitly. 
Eq. ( A. 4 ) is the most general expression for the factorisation at second order in the deriva- 
tive expansion and at finite temperature. 

The temperature dependence of the function B d (£, r, 77) describes dimensional reduction. 
In particular, it obeys the limits 



B d ^,T = 0, V ) =B d {£,Tj) 

B d (t,r >l,r,) = ?±±-*-l-B^ 1 (t,T,) 
v d a — 1 2it 
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(A.13) 
(A. 14) 



The low- and high-temperature limits ( |A.13| ) and ( [A. 14| ) are discussed in the following ap- 
pendices. 

Let us consider the case where £ = 0. It corresponds either to the case of a field- 
independent wave function renormalisation 2 = 1, and/or to the decoupling regime u 3> 1. 
A few properties of Eq. ( |A.4| ) have been discussed in the main text. For £ = the factor 
reduces to 



B d (t,t)) 



d v d -i r 

d — 1 Vd 27T 



2^2' 



1 — 77- 



d + 1 



This corresponds to Eq. ( |4.13| ) discussed in Section [TV]. In addition, we notice that 



fi d (r,r / ^0) = J B,(r) 
^(r-^tj) = 1 

B d (r>l,ri) 



n 



d + 2 

d V d -l T 



d — 1 Vd 2n 



d + 1 



(A.15) 



(A.16) 
(A.17) 

(A.18) 



Eq. ( |A.16| ) corresponds to Eq. (|3.9|), Eq. ( |A.17| ) to the low-temperature limit (|5.15| ), and 
Eq. (|A.18|) to the high-temperature limit ( f4.15|) . 

Similar results are found for the fermionic case, though not discussed explicitly. 



Appendix B. LOW TEMPERATURE LIMIT 

In the low temperature limit r — > 0, the flow ( A.l| ) simplifies to 



%M=fdyy i ? '-\? (1 -f . (B.i) 

J o (2 — 1)2/ + 1 + w 

The remaining integration in Eq. ( B.l|) is solved to give 

4(u,£,ri) = £*(u)B d {Z,ri) (B.2) 

with 

in arbitrary dimensions. For £ = < 1, and hence z > —lu, the integration can be 
peformed analytically. We find 

B d (Z,rj) = B d (Z)-riB i (Z) (B.4) 

where 
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s rf (0 = 2 i ? i(i,i;i + i;e 



(B.5) 



and 



MO = ^ 



^(1,1 + 152 + 1^) 



(B.6) 



For |£| < 1, Eq. ( p.4j ) can be Taylor-expanded in £, leading to the first equation given in 
Eq. (|5TJ). 

For applications, it will be useful to obtain explicit analytical expressions for the functions 
( |B.1|) for fixed dimensions. For d — 4, we find 



1 



1 + UJ Z + UJ 

In 



z-1 0" 1 ) 2 1 + 
-77 



1 + 2uj + z 1 (1 + uj) (uj + z) z + uj 



4(z -l) 2 2 (z-1) 



1 + uj 



(B.7) 



In d = 3 dimensions, we find 

2 2 V / TT^ 



e (u,z,rj) 



{z - 1) - !) 3/2 



arctan 



1 + w 



-77 



1 + 3cj + 2z 
30- l) 2 



1 + Z + UJ 

z-i{z-iy 



arctan 



'z-1 

1 + 



(B.8) 



and the region for z < 1 is obtained through analytical continuation. Finally, for d = 2 we 
find 



1 . z + uj 
In rj 



z-1 1 + UJ 



1 1 



Z + UJ z + UJ 

m 



(B.9) 



2z-l 20 -l) 2 1 + w 

Notice that the functions (|B.7|) , (|B.8|) and ( B.9| ) behave smoothly for z « 1, which follows 
either from Eq. ( p. 14 ) or by an explicit check. 

For z = 1 these expressions reduce to the result Eq. (|5.15|) . 



Appendix C. HIGH TEMPERATURE LIMIT 

The high-temperature limit is reached for r > 1, or T > 2irk. Then, only the n = 
Matsubara mode contributes to the flow. Using Eqs. (|A.7|) , ( |A.10|) and ( P-5|) , we find: 

B d (C, r>l, V ) = -^-^l J-B^fo r = 0, 77) (C.l) 

CL — 1 7J<j Z7T 

For |£| < 1, Eq. ( |C.1| ) can be Taylor-expanded as 

td (c -^.1 \ d ^-1 r » (d-l)f A ry \ 
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Splitting Bd(£, r > 1, 77) as in Eq. (|A.6|) , we have 

u ic ^ i\ d Vd - 1 T 
B d {£, t > 1) = - — — 

a — 1 v d 2it 

In turn, the flow proportional to the anomalous dimension reduces to 

1 d r 



(C.3) 



x 



2d-l v d 2tt 

zji / 1 d— 1 d+1 t\ d-1 771 / 1 d+1 d+3 < 
2-Tl 1 1, — , — ,?J - 7i+T2^1 I — ) — ^ 



(C.4) 
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